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Abstract 


The physical mechanism of fast radio bursts (FRBs) is still unknown. On 2020 April 28, a special radio burst, FRB 
200428, was detected and believed to be associated with the Galactic magnetar SGR 1935+2154. It confirms that 
at least some of the FRBs were generated by magnetars, although the radiation mechanism continues to be debated. 
To this end, we study in detail the multiband afterglows of FRB 200428 described by the synchrotron fireball 
shock model. We find the prediction for the optical and radio afterglows of FRB 200428 is consistent with the 
observations when considering the post-FRB optical and radio upper limits from the literature. We also show that 
the follow up detection of the afterglows from fast radio bursts like—FRB 200428 is possible at the radio band, 
though challenging. Based on our model, one can obtain information about the energy of the fireball, the radiation 
zone, and the nature of the surrounding medium. That may shed light on the physical mechanism of FRBs. 
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1. Introduction 


Fast radio bursts (FRBs) are cosmological radio transient 
sources whose physical origins are still under debated (e.g., 
Cordes & Chatterjee 2019; Petroff et al. 2019; Zhang 2020; 
Xiao et al. 2021). On 2020 April 28, a bright radio burst with 
millisecond duration from the Milky Way magnetar SGR 1935 
+2154 was detected simultaneously by the Canadian Hydrogen 
Intensity Mapping Experiment (CHIME) and Survey for 
Transient Astronomical Radio Emission 2 (STARE2; e.g., 
Bochenek et al. 2020; CHIME/FRB Collaboration et al. 2020). 
With an energy of ~10°° erg, it is considered likely to be an 
extension of the extragalactic FRBs at the lower energy end 
(e.g., Lu et al. 2020). If this radio burst is indeed an FRB (FRB 
200428), it directly confirms that at least some of the FRBs 
were produced by magnetar flares (e.g., Lyubarsky 2014; 
Beloborodov 2017; Margalit et al. 2018; Metzger et al. 2019). 

Due to the more comprehensive data available on repeated 
FRBs, some theoretical models about the FRB generation 
mechanism are proposed. Beloborodov (2017) explains the 
mechanism of repeated FRBs, he suggests that a young 
magnetar star releases energy from successive flares, the flares 
are driven by accelerated ambipolar diffusion in the neutron 
star core, and then power the nebula particles to produce bright 
millisecond bursts. The energy supply of successive flares and 
the collision between different flares give rise to a series of 
repeated FRBs with different intervals. This powerful shock wave 
can also produce bright optical radiation (e.g., Beloborodov 
2020). Yang et al. (2019) studied the brightness and detection 
prospect of “fast optical bursts” (FOBs) associated with FRBs, 


indicating that it is possible to detect the associative FOBs in 
some special inverse Compton scattering processes and with 
telescopes which with high sensitivity. However, Waxman 
(2017) put strict constraints on the nature of the persistent source 
and found that the typical magnetar star wind nebula model is not 
consistent with the predicted results and a strong synchronous 
maser emission mechanism adapted to the GHz band is proposed. 
However, Waxman (2017) found that the typical magnetar wind 
nebulae model did not agree with the predicted results by placing 
strict constraints on the properties of the persistent source, so they 
proposed that the strong synchronous maser emission mechanism 
could be adapted to the GHz band. Metzger et al. (2019) 
suggested that the repeated FRBs might be formed by the 
magnetar-powered synchrotron maser shock model, where the 
central engines released clean ultrarelativistic magnetization 
shock waves spread outward, and then the shock colliding with 
the upstream mildly relativistic magnetized ion-electron shell. 
The shell through the reverse and forward shock wave to 
decelerate, the latter of which produces the observed FRBs by 
synchronous maser mechanism. Similar to GRBs afterglow, the 
forward shock also heats electrons to extremely relativistic 
temperatures, powering (incoherent) synchrotron X-ray /gamma- 
ray emission. Recently, Cooper et al. (2022) have predicted the 
multi-band afterglow of FRB 200428 based on the magnetar- 
powered synchrotron maser shock model. In particular, the 
mechanism of FRB 200428 may be consistent with the 
previously predicted model, but due to its lower luminosity, the 
energy distribution may be slightly different (e.g., Wang et al. 
2020). 
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Another breakthrough in the observation of FRB 200428 
was that the X-ray burst associated with it was simultaneously 
detected, with an energy ratio of the radio burst to the X-ray 
bursts ~10-° (e.g., Mereghetti et al. 2020; Li et al. 2021; 
Ridnaia et al. 2021; Tavani et al. 2021). Incredibly, the far- 
away model has predicted the occurrence of X-ray burst 
associated with the FRB, as well as the low radiation efficiency 
of the FRB (~10-°), which was well verified in the case of 
FRB 200428 (e.g., Margalit et al. 2020). However, recently Wu 
et al. (2020) found that the ejected baryon matter of the 
magenta mainly provided by the crust, is higher than the typical 
mass of a magnetar outer crust. This finding indicates that the 
outer crust of the magnetar predicted in their model cannot 
eject enough baryonic mass. We note that the physical 
mechanism and the rate of the baryonic mass ejection remain 
uncertain. Thus we suggest that further investigations in 
observation and theory are required. The close-in model can 
also account for the low radiative efficiency of FRB 200428 
and its associated X-ray bursts, although they were not 
predicted by the model (e.g., Lu et al. 2020). 

Surprisingly, despite numerous X-ray bursts during the SGR 
1935+-2154 active period, no other FRB event has been detected 
so far, except for FRB 200428 (e.g., Lin et al. 2020; Bailes et al. 
2021; Kirsten et al. 2021). One noticed that the spectrum of the 
X-ray burst associated with FRB 200428 seems to be much harder 
than other bursts, with a peak energy ~85 keV (e.g., Ridnaia et al. 
2021). This may imply that the X-ray burst associated with FRB 
200428 was an unusual burst, and FRBs were produced only in 
this kind of bursts. It thus leads to speculation that this may be the 
underlying reason why other X-ray bursts do not have FRB 
associations. As suggested by Ioka (2020), the usual X-ray bursts 
could come from the fireballs trapped in the closed field lines of 
the magnetar, in analogy to the standard model for soft gamma- 
ray repeaters (SGRs; e.g., Thompson & Duncan 1995, 2001; 
Kaspi & Beloborodov 2017). However, the X-ray burst associated 
with the FRB 200428 may come from the trapped-expanding 
fireball located in the open field line region of the magnetar (e.g., 
Ioka 2020). The observed temperature of such the expanding 
fireball remains constant due to relativistic effects, which is 
consistent with the burst having a larger E, ~ 85 keV than the 
usual bursts. 

As one can see in the physical picture above, the X-ray 
bursts associated with FRB could be accompanied by energetic 
ejecta (the expanding fireballs), while the usual X-ray bursts are 
not. According to the standard fireball model, a pair of shocks 
will be generated when the fireball sweeps the surrounding 
medium, including the reverse shock that propagates through 
the ejecta and the forward shock that propagates through the 
surrounding medium. Similar to the case of gamma-ray bursts 
(GRBs), these shocks would produce multiband afterglows 
(e.g., Yi et al. 2014). The evolution of the afterglows are 
closely related to the energy of the fireball and the nature of the 
surrounding environment (e.g., Yi et al. 2014; Zhang 2014). 
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Therefore, when FRB 200428 repeats again in the future, if its 
multiband afterglows are observed, it will provide new insights 
into the physical mechanism of FRBs. 

In this work, we study in detail the multiband afterglows of 
FRB 200428, and its detectability. It can be seen that there are 
many models that attempt to explain the origin of FRBs. It is 
too early to say which is correct, given the observational 
evidence so far. Therefore, this work does not attempt to 
discuss the framework of any specific FRB model. We refer to 
the afterglow model of GRB and use a memoryless fireball to 
calculate the afterglow of fast radio bursts. Based on such a 
model-independent afterglow model, on the one hand, one can 
predict the timescale and brightness of the afterglow of FRB 
200428 that may occur again in the future according to several 
limited parameters, so as to provide a theoretical basis for 
afterglow observation. On the other hand, future afterglow 
observations can be based on our model to make model- 
independent constraints on the surrounding environment of 
FRB 200428, which is critical for revealing the physical 
mechanism of FRB 200428. 

This work is organized as follows. The standard afterglow 
external shock model is described in Section 2 and the results of 
multifrequency afterglows for FRB 200428 are shown in 
Section 3. The discussion and conclusion are listed in Section 4. 


2. The Model 


Following the previous study of Yi et al. (2014), we applied 
the standard afterglow external shock synchrotron emission 
model of GRBs (e.g., Mészáros & Rees 1997; Sari et al. 1998; 
Gao et al. 2013; Yi et al. 2014) to FRBs. This model describes 
the interaction between the outflow and the ambient medium, 
and it has several free parameters: the total kinetic energy E, the 
number density of the ambient medium nọ, the initial Lorentz 
factor I’, the shock energy equipartition parameters £, and €g 
for electrons and magnetic fields, respectively, and the electron 
injection spectral index p. If the forward and reverse shocks are 
both considered, the equipartition parameters €, and €g and p 
may be different, so there will be nine parameters. 

We mainly consider the forward shock (FS) emission. 
However, an uncertain phenomenon is whether the reverse 
shock (RS) emission appeared, which depends on a magnetiza- 
tion parameter of the outflow o, and the parameter is the ratio 
between the Poynting flux and the matter flux (e.g., Zhang & 
Kobayashi 2005; Mimica et al. 2009; Mizuno et al. 2009). 
Considering FRB 200428 originates from SGR J1935+2154, 
the outflow may be magnetized. As the study of the previous 
(e.g., Komissarov et al. 2009; Granot et al. 2011), the outflow 
is accelerated by a magnetic pressure gradient, and o is 
decreased with the radius, lead to the parameter I’ is increased 
with the radius. Additionally, there is also an obvious magnetic 
dissipation in the phase of the FRB emission. Therefore, the 
value of ø has a lot of uncertainty. If it is already below unity, 
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the RS emission must be expected (e.g., Zhang et al. 2003; 
Zhang & Kobayashi 2005; Yi et al. 2014). 

For the sake of simplicity, we simply considered the standard 
synchrotron emission, which is mainly decided by the bulk 
Lorentz factor [ and the total kinetic energy E, and neglected 
other complication factors. Under this model, the light curves 
evolution are associated with three characteristic frequencies: the 
minimum synchrotron frequency V, the cooling frequency ve, the 
self-absorption frequency v,, and the peak flux of the spectrum 
F; max (€.g., Gao et al. 2013; Yi et al. 2014). The deceleration 
timescale of the thin shell case can be represented as, 


Id +z) 
2cP8/3 ” 


where | = (3E /4rnom,c?)!/ is the Sedov length. 

Although the surface dipolar magnetic field strength of the 
magnetar SGR J1935+2154 is around B, ~ 10G, the radiation 
efficiency of FRBs in radio band is lower, and the isotropic energy 
of FRB may be about 5—6 orders of magnitude less than the total 
energy. According to the observations, the isotropic energy of 
FRB 200428 is about ~10%” erg, therefore we apply the total 
kinetic energy E to be 10% and 10*' erg in this work. 
Additionally, the bulk Lorentz factor of FRB 200428 is adopted 
as a conservative value 50. Even if we set the Lorentz factor to 
100, the brightness of the multi-band afterglows almost has no 
change (e.g., Falcke & Rezzolla 2014; Katz 2014). The 
deceleration time is much shorter than the general GRBs, 
applying the parameters of E = 10“ erg, no = | cm *, T = 50, 
the deceleration time of FRB 200428 is approximately ~0.1 s. 

Based on the standard assumption (e.g., Sari et al. 1998; Wu 
et al. 2003; Gao et al. 2013; Yi et al. 2013, 2014, 2020), we can 
calculate the FRB afterglow emissions. As FRB 200428 is a 
Galactic transient with the luminosity distance of ~10 kpc 
(e.g., Zhong et al. 2020), the FS emission of FRB 200428 can 
be represented at the shock crossing time t,, as, 


(1) 


x 


k 4.1 x 10? ezf -1n T3 Hz, (2) 

=7.5 x 102e Pues ioe E O Hz, (3) 
vf = 4.1 x a Ez}, -1ng Bay Hz, (4) 
Ff maxx = 8.2 x 106%? ono En Dr Jy, (5) 


where €.¢= 0.1, €gf= 0.01 and p = 2.5. The evolution of the 
four parameters during the RS crossing phase (e.g., Mészáros 
& Rees 1997; Sari et al. 1998; Gao et al. 2013; Yi et al. 
2013, 2014), 


t< te vl x 13, vt œ 19, vf ox t 2, FS max X t’, (6) 
and after the RS crosses, 
t > tx: vf œ t, vf, x 177, vf x t3, Ff nax x £9, (7) 


Different from GRBs, due to the FRBs total energy is lower, 
the outflow will rapidly reach the non-relativistic phase, so the 
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dynamics evolution will have different forms. The transition 
time appears when the bulk Lorentz factor [—1=1, here 
T ~ GE/32nnom,c*t?)'/8. The evolution of the four para- 
meters in this phase can be shown as, 
vi ax t5, vf œ t°, vf œ t5, Ff nax x t5. (8) 

The transition time at the non-relativistic phase for FRB 
200428 afterglow is roughly ty ~ 306.1s and 659.5s with 
E=10*° and 10*! erg, respectively. 

If an obvious RS emission exists, the four parameters 
@ and F; max Of the RS at tx can be written as: 


"es = 13 x 103ef 7 1E% -10 T3 Hz, (9) 


r 
Vm VY 


Vi y= 24 x 102e ay TI Eg Hz, (10) 


vi = 4.5 x ita e no^ TSE Hz, (11) 


FT 
v,max, X 


= 2.6 eY 2an DED; Jy. (12) 


Due to the outflow is magnetized, the €g, is normalized to 0.1. 
The evolution of the four parameters of RS before and after 
the crossing time can be shown as, 


t< tg V} oc in, vt vi et, Fn Oct, (13) 


v,max 


and 


t > bxi v x tf T5, vi x tf 3, v’. OC tT 3, F} max X 15. (14) 


2.1. Results 


As the model prediction calculated, the theoretical multi- 
wavelength afterglows of FRB 200428 are shown in Figures 1 
and 2, respectively. Figure | shows the FS afterglow light 
curves of FRB 200428 in the X-ray (2 keV, panel (a)), optical 
(R-band, panel (b)), and radio (1 GHz, panel (c)) bands, 
respectively. Considering different emission efficiencies, we 
take two different kinetic energy values, i.e., E= 10° er 
(black) and 10"! erg (blue). The other parameters are taken as 
the typical values: [= 50, no = 1 cm~°?, Eef=0.1, €g f= 0.01, 
and p = 2.5. Then, we also plotted the sensitivity lines of four 
different detectors in different energy bands, as mentioned in 
our previous study (e.g., Yi et al. 2014). The rad dash line in 
panel (a) is the detector sensitivity line of Swift/XRT, which 
is œf! early on and breaks to ot /? when F,= 
2.0 x 10° ergem™*s"' at t=10°s (e.g., Moretti et al. 
2009; Yi et al. 2014). The green dashed line in panel (a) is 
the sensitivity line of the Insight-HXMT, which scales as at! 
for arbitrarily long exposure times. The red dashed line in panel 
(b) is the sensitivity line of the Large Synoptic Survey 
Telescope (LSST) Array. In the survey model, LSST reaches 
24.5 mag in 30s. The red dashed line in panel (c) is the 
sensitivity line of the Expanded Very Large Array (EVLA), 
which scales as xf !/ for arbitrarily long exposure times. 
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Figure 1. Example forward shock afterglow light curves of FRB 200428. The model parameters: eg s= 0.01, € = 0.1, no = 1, p = 2.5, and n = 50. Three values of 


the energy E = 10°° (black), 10*! (blue) and 10*° (green, only in the optical band) have been adopted. (a) The X-ray light curves at 2 keV. The red dashed line is the 
detector sensitivity line of Swift/XRT, and the green dashed line is the detector sensitivity line of Insight-HXMT. (b) R-band light curves. The red dashed line is the 
detector sensitivity line of LSST. The orange line is the optical upper limits taken from Cooper et al. (2022), and the blue (R-band) and orange (z-band) triangle points 
are the optical upper limits from LCOGT and BOOTES-2 (Lin et al. 2020), respectively. (c) Radio light curves at 1 GHz. The red dashed line is the detector sensitivity 


line of VLA. The green (1.36 GHz) and red (6 GHz) points are the radio upper limits by Effelsberg and VLA (Bailes et al. 2021), respectively. 
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Figure 2. Example reverse shock afterglow light curves of FRB 200428. The model parameters: €g , = 0.16, € = 0.1, n l, 2.5, and 7 = 50. The ener 
p 8 g P B.r e 0 P f gy 


E = 10"' (blue solid line) has been applied. The different detector sensitivity lines are the same as shown in Figure 1. 


As shown in Figure 1 panel (a), the X-ray afterglow of FRB 
200428 is too faint, theoretically neither the Swift/XRT nor 
Insight-HXMT can detect its X-ray afterglow. Considering the 
first X-ray counterpart of FRB 200428 was observed, Insight- 
HXMT implemented a long time observation of SGR J1935 
+2154 since then. Still, no X-ray counterparts from FRB 
200428 were detected, only hundreds of short X-ray bursts 
triggered by Insight-HXMT and a series of other astronomical 
satellites were obtained (e.g., Cai et al. 2022). In the optical R 
band (panel (b)), the peak magnitude is about 16 and 13.6 for 
the energies E = 10°° and 10°! erg, respectively. We can see 
the sensitivity line of LSSA is below the peak magnitude. 
However, due to the very early peak time of ~0.2—0.4 s, which 
lead the LSSA cannot follow up quickly. For optical bands, it is 
still possible to detect the counterpart emission if follow-up 
observations are performed within a few hundred seconds after 
FRB occurs. In the 1 GHz radio band (panel (c)), the peak flux 
density is about 1.46 x 107° and 1.46 x 10°? Jy, the peak time 


is about ~5.7 x 10° and 1.2 x 10° s for two different kinetic 
energies, respectively. This might be caught by EVLA if 
followed up early. Regrettably, the EVLA did not detect this 
source during this period. As reported by Bailes et al. (2021), 
after about 4 x 10* s of FRB 200428 triggered, the radio 
telescope MeerKAT started to point to the source. However, 
they did not detect any signals, due to the flux of the FS 
emission declining rapidly or the diffuse emission around the 
magnetar. Compared with other energy bands, the radio 
afterglows from FRBs are the most promising to be detected, 
mainly because the radio afterglows last a very long time. 
Therefore, if we are lucky enough to detect a bonafide FRB 
like-FRB 200428 again, rapid radio observations on minutes to 
several hours timescales will be the best opportunity to observe 
the afterglows. 

The RS emission for FRB 200428 is also shown in Figure 2. 
Fixing other parameters, we allow the total energy E = 10"! erg 
and £g „= 0.16. In general, the RS afterglow of FRB 200428 is 


Research in Astronomy and Astrophysics, 23:115010 (6pp), 2023 November 


Du et al. 


10 ns 


E(erg) 


i 1 
No» 


logi0(Flux)(Jy) 


1 
w a 


wooo 


Figure 3. Contour of optical peak flux in the E-ng plane. Different colors mean different optical peak flux, and the solid line represents peak flux for 0.1 Jy which is 
the BOOTES-3 upper limit. The black line represents the limitation on the model parameter space imposed by the observation upper limit of BOOTES-3. 


more difficult to be detected by current detectors. Either the 
afterglow emission is too faint, or it peaks too early. As shown 
in Figure 2, like the FS emission, the RS X-ray afterglow 
(panel (a)) is too faint, however the peak time of the optical R- 
band (panel (b)) is so early (t,~0.1 s) for the LSST. 
Additionally, the 1 GHz radio band of FRB 200428 (panel 
(c)) can reach the peak flux early (~33 s) but rapidly declines. 

According to Cooper et al. (2022), who have provided 
optical lightcurves using the method from Margalit et al. 
(2020), and applied the BOOTES-3 upper limits to constrain 
the FRB afterglow significantly. They supposed that the optical 
flux limit scales as Fyjmit X t! and started with early 
observations (the orange line). We also used the same energy 
E=1x 10° erg in our afterglow model, the estimated results 
(the green line in panel (b) of Figure 1) almost reach the optical 
upper limits. Our optical results are also consistent with the 
prediction from Cooper et al. (2022). To better show how 
optical peak flux depends on E and no, we give the contour of 
optical peak flux in the E-ng space in Figure 3. We set 
€g s= 9.01, €e = 0.1, p =2.5 and 7 = 50, but the energy with a 
range of 10*°-10* erg and no with a range of 10 *-10° cm >. 
The reasonable parameter space of FRB 200428 may be 
reflected by the optical observation upper limits from the 
contour of peak flux. As shown in the figure, the black line 
represents the limitation on the model parameter space imposed 
by the observation upper limit of BOOTES-3. It is easy to see 
from the contour that only when the environment around FRB 
200428 is dense enough, it is possible to produce sufficiently 
bright radio afterglows. 


3. Summary and Discussion 


In this work, we applied the standard afterglow external 
shock synchrotron emission model of GRBs to the peculiar 
case of FRB 200428, and calculated its multi-wavelength 
afterglows. As a result, we found that owing to its low energy, 


the broad-band afterglows of FRB 200428 are very faint. Even 
so, current detectors may be able to follow up and detect its 
broad-band afterglows, especially at radio wavelengths. 

The X-ray afterglow emitted from FS is so weak and decays 
rapidly after the peak that it is hard for current detectors to 
catch. This may require a wide-field XRT (such as the Einstein 
Probe or Lobster), and it may have a chance of capturing the 
X-ray afterglow of FRB 200428 in the future, when FRB 
200428 repeats. For the optical afterglow, the current detectors 
can only pick it up theoretically, because it is difficult to follow 
up quickly due to the early peak and the fast decay. The most 
optimistic afterglow of FRB 200428 was found in the radio 
band, although the radio afterglow flux only reached the mJy 
level. Thanks to the relatively late peaking time of the radio 
afterglow (~10°s), if the radio telescopes can move to the 
source position within an hour, there is a good chance of 
catching it. In general, follow up detection of the afterglow 
from the FS of FRB 200428 is plausible, though challenging. 
In contrast, the afterglow from RS is almost impossible to be 
detected by current detectors, either the emission is too weak or 
it peaks too early. Given the difficulty of follow-up observa- 
tion, the strategy of long-term observation when SGR 1935 
+2154 enters the active phase may be adopted. In this way, the 
chance of catching the afterglow is much greater, but much 
more observation time is required. 

Cooper et al. (2022) also estimated the multi-wavelength 
afterglows of FRB 200428, based on the model of Metzger 
et al. (2019). They provided results from LOFAR imaging 
observations of SGR 1935+2154. Due to its low luminosity, 
the predicted multi-band afterglows from FRB 200428 are still 
too faint to be detected, and placing some radio and optical 
upper limits in the afterglow emissions (also see the upper 
limits in our results). Considering the similarity of the early 
afterglow models, our results are very consistent with the 
predictions of Cooper et al. (2022), especially considering the 
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same parameters for this burst. In any case, one can imagine 
that once the afterglow of fast radio bursts like-FRB 200428 is 
observed in the future, based on the early afterglow model, one 
can obtain information about the energy of the fireball, the 
radiation zone and the nature of the surrounding medium. That 
may shed light on the physical mechanism of FRBs. 

Note: Interestingly, since 2022 October 10, SGR J1935 
+2154 is active again (e.g., Ryder et al. 2022) and a radio burst 
was detected from it associated with an X-ray burst (e.g., Dong 
& CHIME/FRB Collaboration 2022; Wang et al. 2022), 
proving another case that magnetars can drive FRBs. In order 
to search for additional FRB-magnetar burst associated cases 
and search for radio pulsations (e.g., Zhu et al. 2020), four 
FAST observations have been done on SGR J1935+2154, but 
no pulses are detected. We also propose four NICER 
observations on SGR J1935+2154 simultaneous with FAST, 
to cover a possible magnetar burst in soft X-ray, and search for 
possible correlated radio/X-ray pulsation, still no pulses are 
detected. On 2022 October 22, a radio burst was detected (e.g., 
Huang et al. 2022) from S-band 40 m Yunnan telescope, CAS. 
Time of S-band burst (MJD): 59873.41789 1687764495, 2020- 
10-21 18:01:45 (UT+8), this arrival time of the radio burst is 
well consistent with an HXMT X-ray outburst (e.g., Li et al. 
2022), the arrival time delay between radio and X-ray can be 
ignored within the error. There is no counterpart detected at this 
stage, and the follow up detection of the multi-band afterglows 
from FRBs is challenging. 
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